Strongly correlated gases of Rydberg-dressed atoms: quantum and classical dynamics 
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We discuss techniques to generate long-range interactions in a gas of groundstate alkali atoms, by 
weakly admixing excited Rydberg states with laser light. This provides a tool to engineer strongly 
correlated phases with reduced decoherence from inelastic collisions and spontaneous emission. As an 
illustration, we discuss the quantum phases of dressed atoms with dipole-dipole interactions confined 
in a harmonic potential, as relevant to experiments. We show that residual spontaneous emission 
from the Rydberg state acts as a heating mechanism, leading to a quantum-classical crossover. 
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There is currently significant interest in the physics 
of dipolar quantum degenerate gases Strong long- 
range dipole-dipole interactions promise the realization 
of novel many-body phases in neutral gases, such as self- 
assembled crystals [2], topological superfluids and quan- 
tum phases with hidden topological order Q . The regime 
of strong dipolar couplings is easily accessible with inter- 
acting electric dipole moments, as realized in particu- 
lar in quantum gases with polar molecules prepared in 
their rovibrational ground state In contrast, quan- 
tum gases of ground state atoms typically interact via the 
much smaller magnetic interactions [Hj]. The question, 
therefore, is to what extent this regime of strong dipo- 
lar interactions can also be realized with present atomic 
gases experiments with alkali atoms. Here we propose 
and investigate a setup where the huge electric dipole 
moments d ~ n 2 [6|-[8| of atomic Rydberg states with 
principal quantum number n are weakly admixed to the 
atomic ground state, thus providing an atomic gas of in- 
teracting effective electric dipoles comparable to the case 
of polar molecules. A central question, and the main dif- 
ference to the molecular case, is decoherence and heating 
mechanisms associated with spontaneous emission from 
Rydberg states, and possible inelastic collisions. Below 
we show that (i) dipolar crystals can be realized with 
Rydberg-dressed atoms confined to two dimensions, with 
negligible collisional losses, and (ii) spontaneous emission 
7 r provides an intrinsic heating mechanism, allowing the 
study of thermalization phenomena with full tunability 
of system parameters. We study the associated crystal 
melting with molecular dynamics calculations. 

The setup we have in mind is illustrated in Fig. [T] we 
propose to weakly couple with laser light the ground- 
state \g) of each atom to a Stark-split Rydberg state 
\r) with large dipole moment cfo, in the kDebye. For 
large enough detuning A from resonance and interpar- 
ticle distances, interactions are of the dipole-dipole type 
oc (Q/A) 4 dl/r 3 , with Q the Rabi frequency. By 
confining the particles to a 2D plane using an optical 
field, the effective in-plane interactions are then 




FIG. 1: (color online) (a) Sketch of the energy levels constitut- 
ing the Stark fan of an alkali metal atom exposed to an electric 
field. \g) and \r) are coupled by a laser with Rabi frequency Q 
and (blue) detuning A. j r is the spontaneous emission from 
|r), with momentum Jcr of the photon recoil, (b) BO poten- 
tials in the (x — ?/)-plane in the dressed picture. Ro is the 
resonant Condon point. The effective interaction potential 
Vjn?( r ) is the higher-energy curve (red). Level crossings oc- 
cur for R n < Ro. Inset: blow-up of V^(r) (continuous line) 
compared to 1/r 3 (dashed line), (c) Sketch of experimental 
setup: Rydberg-dressed atoms are confined to 2D by a strong 
confining laser beam, with dipoles polarized perpendicular to 
the plane. An effective spontaneous emission with rate j e s 
provides for an intrinsic heating mechanism (see text). 



purely repulsive, with negligible collisional losses. This 
opens the way to the study of the many-body phases 
of 2D dipoles in these systems. As an illustration, we 
show the existence of mesoscopic supersolids and crystals 
with Rydberg-dressed atoms under realistic conditions 
of in-plane harmonic confinement, using exact quantum 
Monte- Carlo simulations. Residual spontaneous emis- 
sion 7 e ff ~ (^/A) 2 7 r from the Rydberg state introduces 
an intrinsic heating mechanism, driving the quantum 
phases into the classical regime. We study the quan- 
tum/classical crossover by means of molecular dynam- 
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ics simulations, and show the emergence of a dynamical 
thermalization timescale in these systems. 

Effective interaction potentials with suppressed deco- 
herence from atomic collisions and spontaneous emission 
are obtained as follows. In the presence of a homoge- 
neous electric field (field strength F, oriented along e z ) 
the energy levels of an alkali atom show the well-known 
Stark structure @. We are interested here in high an- 
gular momentum states which show a linear Stark ef- 
fect. The energy of these states is well approximated by 
£nmn 2 ~ (3/2)eaon(ni — n2)F 1 with ri\ and n2 parabolic 
quantum numbers, Oq the Bohr radius and e the elec- 
tron charge. The highest-energy state \r) of a given 
manifold n (n\ = n — 1, = 0) is energetically sep- 
arated from the next (lower lying) state in the mani- 
fold by AE n = e(3/2)aonF, while adjacent manifolds 
remain well separated by an energy AEi T ^> AE n for 
field strengths smaller than Fit ~ 3^? i2gwr 3 e 3 fi 4 (I n gli s- 
Teller limit). In our scheme, Fig. [TJ each atom is ef- 
fectively reduced to a two-level system. Within the 
rotating-wave approximation and with a product ba- 
sis ^\r) i \r) j , \r) i \g) j , \g) i \r) j , I^J^-} the Hamiltonian 
governing the internal dynamics of two atoms is given by 
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with Vij = D(l — 3 cos Oij) |R^ — Rj|~ 
dipole-dipole interaction between two atoms at position 
Ri(Rj), and D = ^Q/(47reo). % is the angle between 
the vector R^ — Rj and the dipole moment aligned par- 
allel to the z-axis. The interaction potential which is 
experienced by the dressed ground state atoms is given 
by the Born-Oppenheimer (BO) energy surface that adi- 
abatically connects to the energy of the product state 
\g)i\g)j as A ^ 00. Here, we focus on the weak driv- 
ing limit Q < A. For two atoms approaching each other 
and blue detuning A > 0, the groundstate BO potential 
will be approximately given by V^(Rij, #;</) c± Duij/Rfj 
for Rij < Ro, with Ro ~ (D/hA) 1 / 3 a resonant Condon 
point with typical values in the hundreds of nm, Fig.[TJb). 
Diabatic crossings with different potential surfaces of the 
same and different n-manifolds leading to collisional two- 
body losses will occur in this region R^ < R , at distances 
of order R n ~ (D/AE n f^ and R IT ~ (D/AE IT ) 1 / 3 , 
respectively, with Rq > R n ^> Rit, reminiscent of blue- 
shielding techniques [2|, [9[ . 

For R^ > Ro, two-body diabatic losses are absent and 
V£?(iJij,0tj) ^ (n/A^Duij/Rlj. In this work we will 
focus on this parameter regime. This has two additional 
advantages: (i) spontaneous emission rates are strongly 
reduced to values 7 e ff and (ii) effective interactions are 
reduced to values compatible with trapping of atoms with 
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FIG. 2: (color online) (a-d) Monte Carlo snapshots of the den- 
sity of particles in all mesoscopic phases for N = 13 dipoles, 
as a function of the effective mass r. (a) superfluid; (b) su- 
persolid; (c) ring-like crystals; (d) classical crystal. 



optical fields (e.g., optical lattice), and confinement to 
low-dimensional geometries. 

Collisional losses for R^ > Ro are linked to popula- 
tion of the attractive part of the dipole-dipole interac- 
tion. Sampling of this attractive part can be suppressed 
by confining atoms to 2D using a tight (optical) trap- 
ping along z, with harmonic frequency huo± = h 2 /ma 2 ± > 
V^(R), with m the atomic mass, and a± the transverse 
harmonic oscillator length. The BO potential is then 
purely repulsive in 3D: there is an energy barrier between 
the long-distance repulsion and the attractive short- 
distance regime. Residual losses are linked to tunneling 
below this energy barrier, and can be computed semiclas- 
sicaly as r co ii = oj at t ex-p[—c(Q 4z Dm/h 2 A 4z a±) 2 / 5 }, with c 
a constant of order unity, and uj att the attempt frequency, 
of order of the average kinetic energy in the gas. For 
Rij > Ro the effective dynamics is purely 2D, with inter- 
actions V 2 ® = (Q/A) 4 D/p 3 = D/p 3 , with p = |p|, and 
p a vector in the (x — ?/)-plane. 

Mesoscopic crystals with Rydberg- dressed atoms: We 
consider a setup where TV bosonic dressed atoms are con- 
fined to a 2D plane by applying a strong transverse trap- 
ping field 0, e.g a ID optical lattice, and are aligned 
perpendicular to the plane [Fig. ffl^c)]. We assume an 
additional in-plane parabolic trap with frequency cj, as 
realized in experiments by a magnetic dipole trap, or 
a single site of a large spacing optical lattice. Defin- 
ing length and energy scales ro = (D/muj 2 ) 1 ^ 5 and 
e = muj 2 ro = D/r^ = (m 3 ^ 6 ^ 2 ) 1 / 5 , respectively, the 
2D Hamiltonian in dimensionless form reads 
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FIG. 3: (color online) Radial density profiles g(r) for the cases 
of Fig. [3] Inset: statistics of computed particle exchanges / 
as a function of the number of particles n' participating to the 
exchange for the cases of finite superfluid fraction r = 1 and 
2.5. In a superfluid, / is finite for all n < TV. The supersolid 
phase with r = 2.5 has finite density modulation (figure) and 
the same /-distribution of r = 1 (inset). 



where r = e/fvuj = (tq/£) = lmD/h 2 £j characterizes 

the strength of the dipole-dipole interactions in the trap 
with i = ^h/muo the harmonic oscillator length. Equa- 
tion (O shows that r plays the role of an effective mass 
representing a control parameter, which can be increased 
by increasing the strength of dipole-dipole interactions, 
or by compressing the trap. 

For small r < 1, we expect the kinetic energy to dom- 
inate, and the cluster to be in a weakly interacting su- 
perfluid (SF) phase, while for r ^> 1 the kinetic energy 
becomes negligible, and the system ground state should 
resemble the classical crystalline (CC) configuration ob- 
tained by minimizing the last two terms of Eq. (|2j). We 
obtain an estimate of the critical r c for the crossover from 
the superfluid to the crystal by noting that for a homo- 



geneous system the supernuid-crystal transition occurs 
at tqm = Dm/h 2 a = 18 ± 4 @, where tqm repre- 
sents the ratio of the dipolar interactions D/a 3 to the 
kinetic energy h 2 /ma 2 with a the mean interparticle dis- 
tance. By rewriting r c = (^qm^/^) 2 / 5 , ano ^ approximat- 
ing a ~ we obtain the prediction r c ~ 3, which is 
essentially N- independent. 

We determine the generic zero-temperature phase di- 
agram for trapped dipolar atoms by means of Quantum 
Monte Carlo simulations based on the continuous-space 
Worm Algorithm 11]. We find that for mesoscopic clus- 
ters with TV < 40 four different quantum phases exist: (i) 
a SF for r <C r c , (ii) a mesoscopic supersolid (MS) for 
r < r c , (iii) a ring-shaped crystal (RC) for r > r c , and 
(iv) a CC for r ^> r c . These phases are distinguished by 
measuring the superfluid fraction p s and the radial den- 
sity profile g(r): the SF and MS phases are superfluid 
with p s = 1, while the RC and CC phases have p s = 0; 
the SF has a flat, featureless g(r), while the MS, RC and 
CC phases have significant density modulations. In the 



crystalline RC phase, particles are arranged in concentric 
rings with a fixed number of particles per ring. For small 
enough temperatures, these rings are free to rotate in- 
dependently of each other. In the crystalline CC phase, 
atoms are arranged in the classical crystal configuration, 
at fixed relative positions. Mesoscopic crystals have been 
also found in excitonic materials in Refs.|l2[. 

In the following we provide example results for TV = 13 
atoms, which we found to display all general features of 
mesoscopic clusters with TV < 40. We find that for TV ^> 
40 the system resembles the homogeneous situation, with 
a sharp crossover between SF and CC phases around r c . 

Panels (a-d) in Fig. |3] are snapshots of the particle 
density for TV = 13 and 1 < r < 20 in the four meso- 
scopic phases described above. These results correspond 
to a low enough temperature T <C e, that they can be 
regarded as groundstate estimates. Panel (a) shows a 
SF for r = 1 with a flat density profile (see also g(r) 
in Fig. [3]), where the various particle probability clouds 
overlap. This overlap is directly connected to superfluid- 
ity [13|, and consistently we here measure p s = 1. The 
emergence of a MS phase for r = 2.5 < r c is signaled 
by a distinguishable density modulation in Fig. QJb) (see 
also Fig. [3j), combined to a measured p s = 1. We find 
that the superfluid properties are completely unaltered 
by the increased strength of interactions with respect to 
the case r = 1. This is quantified in the inset of Fig. [3] by 
measuring the statistics of exchange cycles / in a many- 
particle path [13]: in a superfluid, the probability / that 
n' particles exchange is finite for all n' < TV, while in a 
crystal is approximately zero. In the inset of Fig. [3] it 
is shown that / is finite and equal for r = 1 and 2.5, 
computed at corresponding low temperatures (i.e., same 
fractions of e). 

For r > r c [panels (c-d)] we observe crystallization of 
the atomic cloud, with p s = and finite density modu- 
lations, consistently with the estimates above. In partic- 
ular, for r = 5 we obtain a RC phase, with concentric, 
independent rings while for r = 20 particles are ar- 
ranged in the CC configuration for TV = 13, with 4 par- 
ticles at the center, and 9 outside. This classical crystal, 
characterized by large peaks in g(r) in Fig. El is found to 
be the groundstate configuration for all r > 20. 

Effective heating: The residual effective single-particle 
spontaneous emission rate is 7 e ff ~ (Q/A) 2 j r . In this 
work, we focus on the case where after the spontaneous 
emission the atom is found in and then instanta- 
neously re-dressed, which is the relevant situation for 
n < 30 [IH. Spontaneous emission results in the release 
of a photon, with momentum hk and recoil energy Er^ 
in the tens of kHz. For dipoles parallel to F, emission 
will be preferentially along z, and Er will be absorbed by 
the confining lattice potential with Huj± > Er. However, 
any finite component of the momentum hk xy of the emit- 
ted photon in the (x — ?/)-plane will result in a random 
"kick" of the atoms, which will ultimately translate into 
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an effective intrinsic heating rate for the many-body sys- 
tem. We conservatively estimate the crystalline phases 
(r ^> 1) to be observable up to an average "melting" time 
t m oc Tyij '(^QfiEn), with Tm the (classical) melting tem- 
perature in the trap. We computed Tm for all N < 30 
by means of molecular dynamics simulations for the 2D 
system in the classical limit (r — >> oo) [16]. By monitor- 
ing the behavior of the specific heat and density-density 
correlations, we determine the value Tm = 0.4(2)e. This 
can result in characteristic lifetimes of order of several 
milliseconds. For example, for 87 Rb atoms coupled to 
the Rydberg-state (n = 20, ri\ = 19, ri2 = 0) with 
do ~ 1.45kDebye, 7 r ~ 100kHz, Er c± 25kHz, via a laser 
with n/2n = 100MHz, A/2tt = 1GHz and a DC field 
F = 25kV/m, for uj/2tt = 5kHz we obtain tq ~ 680nm, 
r c± 20, e c± 4^, and t m > 5ms. 

We further investigated the melting dynamics. Spon- 
taneous emission was simulated by applying random in- 
plane kicks to each particle with momentum hk xy = 
(2mEn xy ) 1 ^ 2 , with fixed Er xv < Er, at an average rate 
7 e ff. Panels (a) and (b) in Fig. [4] show our results for 
ERxy — 0.5e and 2.5 e as a function of time £, respec- 
tively, for A/" = 13 particles [see Fig. |3]. Thin lines are 
single classical-dynamics trajectories, while thick dashed 
lines are averages over many trajectories. We find that 
the average energy E tot increases linearly with time t, as 
expected. For the smaller Er xv of panel (a), we find that 
the crystal melts roughly at the expected time t m (cor- 
responding to several spontaneous emission processes), 
while in panel (b) the observed melting time t eq is signif- 
icantly longer than t m (roughly the time of one sponta- 
neous emission). This prolonged stability is a dynamical 
process: classical melting can occur only after the system 
has reached thermal equilibrium. The latter corresponds 
to the time t eq at which the average potential energy E pot 
and kinetic energy E^ n become equal and equipartition 
of energy applies. t eq , and thus the crystal lifetime, can 
be (much) longer than t m . 

Since interparticle distances can be of the order of a 
/im or more, the spatial structure of the crystalline phases 




FIG. 4: (color online) Molecular dynamics (MD) simulation 
of crystal heating. Black, red and blue lines are total en- 
ergy .Etot, kinetic i^kin and potential E pot energies vs. time 
t, respectively. Thin continuous lines: single MD trajectories. 
Thick dashed lines: averages over many MD trajectories. t eq 
and t m are the equilibration and melting time, respectively. 



can be imaged using, e.g., tightly focused beams. 
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